Quantum dynamics of a dc-SQUID coupled to an asymmetric Cooper pair transistor 
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We present a theoretical analysis of the quantum dynamics of a superconducting circuit based 
on a highly asymmetric Cooper pair transistor (ACPT) in parallel to a dc-SQUID. Starting from 
the full Hamiltonian we show that the circuit can be modeled as a charge qubit (ACPT) coupled to 
an anharmonic oscillator (dc-SQUID). Depending on the anharmonicity of the SQUID, the Hamil- 
tonian can be reduced either to one that describes two coupled qubits or to the Jaynes-Cummings 
Hamiltonian. Here the dc-SQUID can be viewed as a tunable micron-size resonator. The coupling 
term, which is a combination of a capacitive and a Josephson coupling between the two qubits, can 
be tuned from the very strong- to the zero-coupling regimes. It describes very precisely the tunable 
coupling strength measured in this circuit [l] and explains the 'quantronium' as well as the adiabatic 
quantum transfer read-out. 



I. INTRODUCTION 

Two quantum systems with discrete energy levels cou- 
pled to each other form an elementary block, useful for 
the study of fundamental phenomena and effects in quan- 
tum physics, especially in the context of quantum in- 
formation. The interaction between the two quantum 
systems is essential to implement important concepts in 
this field such as entanglement, quantum gate operations, 
and quantum information transfer. As to the theoretical 
description of interacting quantum systems, two prob- 
lems have been extensively studied in particular: a two- 
level system (or qubit) coupled to a harmonic oscillator 
and two coupled qubits. The former was used to de- 
scribe among others the quantum electro-dynamics asso- 
ciated with atoms in a cavity [2], trapped ions coupled 
to their vibrations and more generally interaction be- 
tween matter and photons pi. The latter was developed 
to describe entangled photons, trapped ions, and two- 
qubit quantum gate operations [5 . 

These studies were experimentally realized in the fields 
of quantum optics and atomic physics and more recently 
extended to include quantum solid state devices. In par- 
ticular, during the last decade superconducting circuits 
have demonstrated their potential in connection with 
quantum experiments (6H15) . They now appear as exper- 
imental model systems for studying fundamental quan- 
tum physics and basic blocks for quantum information. 

In this paper we consider a superconducting circuit 
composed of two well-known elements coupled to each 
other: an inductive dc-SQUID and a Cooper pair tran- 
sistor. This circuit constitutes an elementary building 
block that can be operated in various parameter regimes 
characterized by different types of quantum dynamics, 
as has been shown experimentally in the past. As we 
will detail below, this is possible due to the fact that 
three strongly coupled quantum variables determine the 
dynamics of this coupled circuit. 

For instance, when the current-biased dc-SQUID is 
non-inductive and classical and the transistor is sym- 
metric, one recovers the Quantronium ^6]. When the 



SQUID is operated in the nonlinear regime, the result- 
ing system consists of a charge qubit coupled to an an- 
harmonic oscillator; this system has been shown to al- 
low for non-destructive quantum measurements [I7|. We 
recently operated the SQUID in the quantum few-level 
limit [1^, demonstrating a very strong tunable coupling 
between two different types of qubit: a phase qubit and 
a charge qubit. 

The experimental performance of this circuit was lim- 
ited by uncontrolled decoherence sources. However its 
integration (three quantum variables strongly coupled 
on a micrometer scale), its tunability, its various opti- 
mal points make this circuit attractive once decoherence 
sources will be overcome upon technological improve- 
ments. 

In this paper we present a rigorous theoretical anal- 
ysis of the circuit in the parameter range of our exper- 
iments [1 . The full Hamiltonian of the coupled circuit 
is presented, describing a two- level system (Cooper pair 
transistor) coupled to two anharmonic oscillators (dc- 
SQUID). In the relevant parameter range, the dc-SQUID 
dynamics can be reduced from two-dimensional to one- 
dimensional. Consequently the dynamics of the circuit is 
that of a qubit coupled to a single anharmonic oscillator. 
Depending on the anharmonicity, different regimes can 
be studied in this unique circuit. When the anharmonic- 
ity is neglected, we recover the physics of a qubit cou- 
pled to a harmonic oscillator. Its quantum dynamics can 
then be described by the well-known Jaynes-Cummings 
Hamiltonian. Although this limit can also be achieved 
with a qubit coupled to a high-Q coplanar wave guide 
cavity [18] we wish to emphasize that the use of a 
dc-SQUID is of interest as it constitutes a micron-size 
resonator and it is tunable. When the oscillator is con- 
sidered anharmonic, its interaction with a qubit gives rise 
to very complex dynamics which has been very little ex- 
plored. If only the two lowest levels of the oscillator are 
considered, it can be reduced to a two-level system. The 
coupled circuit then describes two interacting qubits. In 
addition to the possibility to study different dynamical 
regimes depending on the anharmonicity of the resonator, 
the coupling between the SQUID and the transistor is 
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fully tunable. As a result the system can be operated at 
zero coupling, as well as in the weak and in the strong 
coupling limits. 

In section II, after a description of the circuit un- 
der consideration, we construct the Lagrangian and the 
Hamiltonian using Devoret's circuit theory [19]. Section 
III is devoted to the two-dimensional dynamics charac- 
teristic for an inductive SQUID and its reduction to one- 
dimensional dynamics in the relevant parameter range. 
In section IV we discuss the quantum dynamics of the 
asymmetric Cooper pair transistor, especially at its two 
optimal points where it is insensitive to noise fluctua- 
tions. In section V the terms describing the coupling 
between the dc-SQUID and the transistor are discussed. 
In section VI, the full quantum dynamics of the coupled 
circuit is presented. There we also discuss the two pos- 
sible quantum measurements of the charge qubit that 
can be performed by the dc-SQUID. The last section dis- 
cusses the tunable coupling strength of the circuit and 
its comparison with the experiments. 



II. COUPLED CIRCUIT DYNAMICS 

A. Circuit description 

A schematic electronic representation of the circuit 
studied theoretically hereafter is presented in Fig. [ijb). 
In this circuit two different elements are coupled that 
correspond to two basic blocks for typical superconduct- 
ing quantum devices. The first element is a dc-SQUID; 
it corresponds to a loop which contains two identical 
Josephson junctions (JJ), each with a critical current 
Jo and a capacitance Cq. The total inductance Ls of 
dc-SQUID is the sum of three inductances Li, L2 and 
I/3 associated with the different parts of the SQUID 
loop. The second element of the circuit is an asymmet- 
ric Cooper Pair transistor (ACPT) which is placed in 
parallel with the dc-SQUID. The ACPT consists of a su- 
perconducting island connected to the dc-SQUID by two 
different Josephson junctions. We denote by if and /J 
the critical currents and Cf and the capacitances 
of these junctions. The asymmetry of the transistor is 
characterized by the Josephson asymmetry parameter 
/i = {I2 — Ii)/{l2 + ^1^) and the capacitance asymme- 
try parameter A = {Cf - Cj)/{Cf + ). The ACPT 
is also coupled to a gate-voltage; this is modeled theo- 
retically by an infinite capacitance Cp with a charge Qi) 
with Cp^Qb 00 so the ratio Qb/C^ = Vg. The cir- 
cuit is current-biased, modeled similarly with the help 
of an infinite inductance L5 ^ 00 threaded by a flux 
$5 ^ 00 so that the ratio ^5/^5 = h remains constant. 
The properties of the circuit depend on four, experimen- 
tally tunable parameters Vg^ lb and the fluxes and 
$T threading the dc-SQUID loop and the other loop of 
the circuit, respectively. As we will see, these parameters 
allow to control and change the dynamics of the coupled 
circuit. This circuit, realized and studied experimentally 



(a) ACPT dc-SQUID 




FIG. 1. The coupled circuit studied experimentally in Ref. [T]. 
(a) On the right, a SEM picture shows the two elements of 
the circuit, i.e., an Asymmetric Cooper Pair Transistor (red 
frame) connected in parallel to a dc-SQUID (blue frame). A 
zoom of the transistor is shown on the left of the SEM picture. 
The transistor asymmetry stems from the difference between 
the areas of the two Josephson junctions. The quantum dy- 
namics of the circuit can be modified by the bias current 
the flux $s in the dc-SQUID, the flux $t in the other loop 
and the gate- voltage Vg. (b) Schematic representation of the 
circuit. The four Josephson junctions are represented by a 
pure Josephson element in parallel with a capacitance. The 
voltage bias of the transistor gate and the current bias of the 
circuit are schematically indicated by an infinite capacitance 
Cp and an infinite inductance Lp respectively. The spanning 
tree drawn in red connects the ground node to the six active 
nodes cpi, (^2, V^, 0, 7 and ^. 

by A. Fay et al. pLj, is shown in the SEM view in Fig.jlja). 



B. Devoret's circuit theory 

The relevant degrees of freedom of a superconducting 
circuit and their dynamics can be determined using the 
concept of node phases introduced by M. Devoret [19]. 
We distinguish two different kinds of nodes in the cir- 
cuit. We first choose a ground node to which the zero 
of phase is associated. Note that this choice corresponds 
to a choice of gauge and is therefore arbitrary; although 
this choice affects the detailed form of the Hamiltonian 
it does not affect the resulting dynamics of the circuit. 
The other nodes of the circuit are called the active nodes, 
each described by an active phase. Six active phases are 
present in the circuit considered here; they are denoted 
by ipi, (p2^ V^, 0, 7 and ^. 

Let us now introduce the so-called spanning tree [19]. 
Starting from the ground node we draw the branches of 
the spanning tree reaching each active node via a unique 
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FIG. 2. Three dipoles in a loop threaded by the flux The 
spanning tree in red connects the ground node to the two 
active nodes (/)a and (^b- The superconducting phases Ai, A2 
and A3 across the three dipoles respectively are functions of 
4>A, <pB and 



age V is related to the phase difference A by V = ^qA. 
The current Ig through the capacitance is then given 
by Ic = (/>oC'A, where <po = ^o/{27t). The current II 
through the inductance reads II = (/>oA/L. Finally, the 
current through the Josephson element is given by the 
Josephson relation Ij = Iq sin(A), where Iq is the Joseph- 
son critical current |21j. Using these expressions for the 
current, we can write the six current conservation laws 
for each active node in the studied circuit (the sum of the 
currents flowing into each node is zero). The six conser- 
vation laws yield six equations for the dynamics of the 
node phases (cf. appendix [A|), whose solution yields the 
dynamics of the entire circuit. 



D. The Lagrangian 



path. In the case of the coupled circuit, the spanning 
tree (drawn in red in Fig. [ijb)), connects the ground 
node (bottom node) to the six active nodes. 

The superconducting phase difference across a dipole 
in a superconducting circuit can be written as a function 
of the node phases with the help of two rules. We will il- 
lustrate these rules with the help of the circuit presented 
in Fig. [2] as an example. Here, three dipoles are placed 
in a loop threaded by the flux ^. In this case, there are 
three nodes; the spanning tree (drawn in red) connects 
the ground node to the two active nodes with phases (j)A 
and respectively. As a first rule, the superconducting 
phase difference across a dipole located on the spanning 
tree is given by the difference of the phases of the nodes 
linked to the dipole. Hence, the superconducting phases 
Ai and A3 of the dipoles 1 and 3, respectively, are given 
by Ai = (j)A and A2 = 0b- As a second rule, in the case 
of a dipole which is not located on the spanning tree, 
we first define a minimal loop which contains the previ- 
ous dipole and the other dipoles located on the spanning 
tree. The superconducting phase difference across the 
dipole is calculated, using the quantization of the phase 
for the minimal loop [20 . Let us apply this rule to de- 
termine the superconducting phase difference A2 across 
dipole 2. The minimal loop corresponds simply to the 
loop of the circuit. The phase quantization in this loop 
gives A2 = 0B — 0A — 27r<l>/$o, with $0 the supercon- 
ducting flux quantum. With the help of the two previous 
rules, we have analyzed the coupled circuit and expressed 
the superconducting phase difference across each dipole 
as a function of the six active phases. 



C. Current conservation and system dynamics 

In a superconducting circuit, we find generally three 
different dipoles: a capacitance C, an inductance L and a 
Josephson element J. The current through each of these 
dipoles can be expressed as a function of the supercon- 
ducting phase difference A across the dipole. The volt- 



The Lagrangian C of the circuit depends on the six 
node phases and their time derivatives. Let us define 
the vector x = {cpi^ (f2-, V^, ^, 7, formed by the six node 
phases of the circuit. The six Euler-Lagrange equations 
are defined by [22] : 



dx 



0. 



(1) 



The Lagrangian of the circuit has to be constructed in 
such a way that the Euler-Lagrange equations are equiva- 
lent to the current conservation equations ( Appendix [A]). 
We take the Lagrangian to be of the following form: 



,^3 {^x , X ) 

with the kinetic energy 



r(f) - v{x), 



(2) 



+ \cIi^' + \c,{i^-of + \cpe'} (3) 

and the potential energy 

V{x) = -0o{4 cos{(pi) + /o cos((/;2) 

+ /J cos(^ - - (^t) + cos(V^)} 



2 Li 2 L2 

2 Ls ^( 



(4) 



The kinetic energy corresponds to the energy stored in 
the capacitances of the circuit, whereas the potential 
energy is composed of the Josephson energies (cosine 
terms) and the energy stored in the inductances of the 
circuit. From now on, we assume the SQUID inductances 
Li = L2 and introduce the inductance asymmetry r] de- 
fined by L3 = rjLs [23 . Although the Lagrangian de- 
pends on six variables, the effective low-energy behavior 
of the circuit is determined by three variables only, as we 
will see below. 
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Let us first consider the phase variable 7. Its dy- 
namics is that of a harmonic oscillator of frequency 



Lagrangian can then be expressed as: 



node and Ls = r]{l — r])Ls/{l + r]). Using the circuit 
parameters of Ref. [1] (see Appendix [B| , where is es- 
timated to be smaller than 0.1 fF, the frequency uj^ is 
estimated to be larger than 1 THz, i.e., larger than all 
the other frequencies of the circuit. Next consider the 
phase variable ^. Again using the circuit parameters of 
Ref. [1] , we find that the characteristic inductive currents 
^0/^1,2 are on the order of 1 /iA, much larger than the 
Josephson critical current /J of about 1 nA. In other 
words, the SQUID inductance is much smaller than the 
Josephson inductance ^ Therefore, in a first ap- 

proximation, the dynamics of ^ is also of the harmonic 
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1/JLsC^. Here is the capacitance of the 7 T{x,y,^l),e) = (^^<j Cqx'^ + Coy^ ^ {x - rjy - ^j) 
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oscillator type with a frequency uj^ = y l/L^-Cj. We es- 
timate uj^ to be around 640 GHz. This frequency is much 
smaller than cj^, but still quite high compared to the fre- 
quencies characterizing the dynamics of the variables (pi^ 
ip2i and (around 10 GHz, see below). This implies 
that we can use the adiabatic approximation to elimi- 
nate the fast variables 7 and ^ and obtain an effective 
Lagrangian for the slow variables 0i,02,^, and 6. 

In order to implement the adiabatic approximation, 
we write 7 = 70 + ^^7 and similarly = (^q -\- S(_. Here 70 
and ^0 correspond the minima of the harmonic potentials 
confining the motion of these variables, 

-fo{x,y) =x^r]y^ ^(1 -^) + 7^^5^p(1 - ^^). (5) 
2 400 

^o{x,y) =x-r]y- ^(1 -r?) + -^Lslp{l - r]f 
2 400 

where x = {cpi -\- (p2)/2 and y = {(p2 — ^i)/2. Note that 
70 and ^0 depend on x and y^ hence they acquire slow 
dynamics through the phases cpi and (p2. The dynamics 
of the deviations J7 and 5^ is fast, determined by the fre- 
quencies 00 J and cj^, respectively. Substituting the above 
decomposition for 7 and ^ in the Lagrangian ([3|, and 
Q, and averaging over the fast variables (^7 and 5^ we 
find the effective low-frequency potential energy 

V{x,y,il;) 

2Ej (- cos(x) cos(?/) - s{r]y + x) + h{y - ys)'^) 

-EJ^ cos(V^) - cos (V^ - ^o(^, y) + (/>t) . (7) 

Here we defined the reduced parameters b = (j)o/{LIo)^ 
s = h/{2Io) and ys = 7r/{^s/^o)' Furtermore, Ej = 
(polo, EJ^ = (po^i and EJ2 = ^/^o^J are the different 
Josephson energies of the circuit. Note that the fast os- 
cillations of ^ lead to a renormalization of /J; the value 
it takes in the effective low-frequency potential ([7| will 
therefore be smaller than the bare value appearing in 
Q. Assuming that the bias current and the flux 
are constant, we deduce from Eqs. ([5| and ([6| that 
7o = X -\- rjy and ^0 = x — rjy. The kinetic part of the 



+ lcl^' + \c,{i,-e)' + \cpe'y (8) 

The four variables of the Lagrangian can be separated 
in three groups. Indeed, the dynamics of x and y corre- 
sponds to that of the dc-SQUID, whereas the dynamics of 
t/j is associated with that of the ACPT. The variable is 
used to model the effect of the gate voltage (cf. Sec. II F ). 
Note that the last term of the potential ([7| and the third 
term of the kinetic term ([8| couple the variables of the 
dc-SQUID and those of the ACPT together, and there- 
fore are responsible of the coupling between these two 
elements. 



E. Choice of variables for the dc-SQUID 

The Lagrangian of the circuit is a function of the vari- 
ables X and y associated with the SQUID. We want 
to change these variables to more appropriate ones 
which will be used bel ow t o describe the dynamics of 
the dc-SQUID (cf. Sec. HI). Let us introduce the two- 
dimensional potential of the dc-SQUID Vs{x,y) which 
is the contribution to the potential y(x,y,'0), Eq. ([7|), 
depending only on the variables x and y. It reads: 



Vs{x, y) = 2Ej (- cos(x) cos(^) - s{r]y + x) + b{y - ysf) 

(9) 

We stress here that this potential is identical to the one 
of a dc-SQUID alone, as studied by J. Claudon et al. 
[28] . The dynamics of the dc-SQUID is similar to that 
of a fictitious particle of mass ~ (/>o2Co which evolves in 
the potential Vs{x,y). This potential undulates due to 
the product of cosine terms and contains wells that are 
separated by saddle points (see Fig[3|; Vs{x,y) is modu- 
lated by the bias current /5 and the flux ^s- We consider 
now the case that the particle is trapped in one of these 
wells associated with a given local minimum (xo, ^o)- Let 
us introduce the displacement variables around (^0,^0) 
defined hy X = x — xo and Y = y — yo. We assume that 
the particle's motion does not extend far from (xo,^o)- 
Then, we can replace the potential Vs{x,y) by its third 
order expansion around {X = 0, F = 0). This expan- 
sion contains a cross-term in XY which disappears by 
performing a rotation of the (X,Y) plane (Fig. |3jb)) 
by the angle 0, where is given by tan(2^)/2 = 
dlM{xo, yo)/{dl^U{xo, yo) - d^M{xo, yo))' The 



new 



variables of the dc-SQUID Xy and Y± associated with 
the rotated plane are defined by 



^11 



cos{0) sm{0) 
— sin(^) cos(^) 



(10) 



and correspond to the position of the particle along the 
longitudinal and the transverse direction, defined by the 
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y 

FIG. 3. Potential Vs{x,y) of the dc-SQUID with bias pa- 
rameters s = 0.55, Vb = (SQUID parameters b = 1.28, 
rj = 0.29). The point (xo,yo) is the local minimum of one of 
the potential wells. 



two- dimensionality of the SQUID potential. When the 
inductance Ls is zero, the dynamics of the SQUID is de- 
scribed by only one variable, since y = ys- In that 
case, the prefactors ay and a± are equal to one. The two 
last terms of the potential contain one variable of the 
SQUID and one of the transistor. They thus couple the 
SQUID and the transistor. Note that the coupling terms 
of the second order and beyond have been neglected in 
the potential. The kinetic term of the Lagrangian can be 
rewritten as, using the variables Xy and F^, 

T(X||,r^,^,^) = 

-2C^# + (Cp + C,)^2|, (13) 

where Cy: is the total capacitor of the transistor, defined 
by = Cf -\-C2 -\-Cg. The final expression for the total 
Lagrangian is obtained from jC{X\\ , y^, ?/^, ^, Xy , Y^, V^) = 
T(X|| , Fl , V^, ^) - F(X|| , , ?/;) ; it will be used in the next 
section to establish the Hamiltonian of the circuit. 



F. The Hamiltonian 



minimal and the maximal curvature of the potential, re- 
spectively. The third-order expansion of the SQUID po- 
tential now takes the form 



ii-'^ii 



k±Yl- 



a_LY_ 



(11) 



where the prefactors /cy, ay, k_\_, /3a and can be 
calculated numerically. The dynamics of the fictitious 
particle of the SQUID in this potential will be analyzed in 
Sec. |III[ The full potential appearing in the Lagrangian, 
Eq. JtT, is given by: 



y(Xy,n,V^) = 



-fcyX|+ayXff 



-£;Ji cos(V^) - eJ^ cos(Vi - S) 

-Ej^a^Xw sin(^ - S) ^ Ej^a^Y^ sin(^ - S), (12) 

where S = ^(xo,^o) — is the classical phase dif- 
ference across the transistor. The prefactors ay = 
cos(^) — 7^sin(^) and a± = sin(^) + r]cos{0) reflect the 



The Hamiltonian of the coupled circuit is a func- 
tion of the variables (Xy , Y^, ?/^, ^) and the conjugate 
momenta {—hP\\,—hP±^—hn^—hnQ). These momenta 
are related to the velocities (Xy,!^,'^,^) by the well- 
known expressions —hP\\ = dC/dX\\^ —hP±_ = dC/dY^^ 
— hn = dC/d^ and —Huq = dC/dO. The analytical 
expressions for the conjugate momentum variables are 
given in appendix [C| The conjugate variables generate 
the charges Py, n and uq with unit [-2e]. We stress that 
these charges have a clear physical meaning. Indeed, Py 
corresponds to the number of Cooper pairs stored in the 
two capacitors Co; n is the number of Cooper pairs on 
the island. In Eq. (C4), the charge —2enQ is equal to 
the bias charge Q^. Performing the limiting procedure 
Cp^Qb oo, keeping their ratio constant Qb/Cp = Vg, 
we see that the velocity is constant and defined by 
= Vg/(j)o. The expression for the Hamiltonian is deter- 
mined by the Legendre transformation [22]: 



-hPuXu - hP±_Y^ - hnijj - C. 



(14) 



All the velocities which appear in the Lagrangian have to 
be replaced by the conjugate variables, inverting the set 
of equations given in appendix IC] The full Hamiltonian 



of the circuit then takes the form: 



(2ef 
2C± 

C||± ' 
{2ef 

C'nW 



2e 



£;Ji cos(?/^) - cos(^ - 5) 



^)-^J,a||X||sin(^-J) 

CaV, 



n - 



2e 



■Ej^a^ sin(^-(5), (15) 



where the analytical expressions of the capacitances Cy, 
C±, Ct, Cn\\, Cn± and C\\j_ are given in Tab. Ill Apply- 
ing the standard canonical quantization rules, the classi- 
cal variables have been replaced by their corresponding 
quantum operators. For more clarity, the conjugate pairs 
appear in the Hamiltonian in different colors. They sat- 
isfy the following commutation operations 



(16) 




The properties of the Hamiltonian (15) are not trivial. 



It describes the quantum dynamics of three sub-systems: 
the longitudinal and transverse phase oscillations within 
the dc-SQUID and the charge dynamics of the ACPT. 
Moreover the Hamiltonian describes the dominant cou- 
pling between the different quantum sub-systems. Very 
complex dynamics can appear in this full circuit. In this 
paper we mainly concentrate on the dynamics of the lon- 
gitudinal SQUID phase mode and the charge dynamics 
of the ACPT and as well as on their coupling. The next 
section is dedicated to the study of the dc-SQUID Hamil- 
tonian. We will deduce the simplified Hamiltonian for 
the longitudinal phase mode and justify why transverse 
phase mode can be neglected in this study. 



III. DC-SQUID 



The dc-SQUID potential has already been discussed in 



Sec. |IIE[ where we introduced the change of the variables 
X and y to the variables Xy and Y_\_. In this section, we 
first analyze the properties of the dc-SQUID potential in 
more detail. Then, we study the quantum dynamics of 
the dc-SQUID which is equivalent to that of a fictitious 
particle trapped in one of the wells of the potential. We 
will see under which conditions the dc-SQUID behaves 
as a phase qubit. In this Section, the coupling between 




O /O 
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FIG. 4. Critical current Ic[f] for well family [/] as a function 
of the flux ^s- 



the SQUID and the ACPT will be ignored such that we 
can consider the dc-SQUID as an independent element. 



A. dc-SQUID potential 

For an appropriate choice of bias parameters, i.e., the 
bias current and the flux $5, the SQUID potential 
contains wells that can be regrouped in families [/] [24] , 
The index / for the family [/] is related to the number / 
of flux quanta trapped in the dc-SQUID loop. The wells 
of the same family are located along the direction x, peri- 
odically spaced by a distance 27r, and separated by saddle 
points. Since these wells have exactly the same geometry, 
the physical properties of the SQUID are independent of 
the particular well in which the fictitious particle is local- 
ized. Wells of the family [f] exist only if the bias current 
satisfies the relation /"[/] < h < I^[f]^ where I^[f] 
and /"[/] are the positive and negative critical current 
of the family [/], respectively. When = ^tif]^ 
local minima of the wells of the family [/] and their clos- 
est saddle points coincide. The critical currents depend 
strongly on the flux $5 as shown in Fig. [4] for the ex- 
perimental parameters of the circuit studied in Ref. [1]. 
Here, for a given flux, the absolute values of the critical 
currents I^[f] and /"[/] are different. This difference 
originates from the finite inductance asymmetry of the 
dc-SQUID {f] = 0.28) and disappears when 77 = 0. For 
almost any value of the flux, the potential is character- 
ized by a unique family of wells, except in the region close 
to ±^0/2 where two families can coexist. This specific 
region has been investigated recently in order to study 
the double escape path of the particle, as well as to make 
the SQUID insensitive in first order to current fluctua- 
tions [H]. In the following we will discuss the dynamics 
of the dc-SQUID for the case where the fictitious particle 
is trapped in one of the wells of the potential. 
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capacitance labels 






Ct 


Cn\\ 




^11^ 


exact expressions 


















2Co + Cj(l + r72) 






(a_LQ!|| C_r) 


approximated expressions 






cT + 


2Co(Cj+Cf ) 






numerical values 


455.6 fF 


455.1 fF 


2.90 fF 


653.0 fF 


2.252 pF 


1.17 nF 



TABLE I. 
of G 



Analytical expressions for the capacitances present in Hamiltonian ( 15 ) and the approximated expressions in the 
T p. r^T j^g^^g ^gg^ definitions 



2C7o(Cs + Cr(\ + ry')), Cs = Cl + Cl^ + Cg and 



limit 

= C2 (C^ + Cg)/2CQ. The numerical values have been calculated for the parameters of the circuit studied in Ref. [l] 
Co = 227 fF, Cr = 2.0 fF, = 0.9 fF, Cg = 29 aF, and for a zero escape angle (6' = 0). 



B. Hamiltonian of the dc-SQUID 



as 



In the full Hamiltonian (15) of the coupled circuit, we 



isolate the terms which only contain the operators X\\ 

and and their conjugate momenta P\\ and P^. These 
terms constitute the Hamiltonian of the dc-SQUID which 
takes the form 



Hs = ElPl + -k\\Xl+a\\Xl 



+PaYlX\\ + PkY^Xl + V^PyPi, (17) 



with the charging energies E^^ — (2e)^/(2C||) and = 
(2e)^/(2Cx). The first three terms correspond to the 
Hamiltonian of a fictitious particle of mass m = (i)\C\\ 
which is trapped in an anharmonic potential along the 
longitudinal direction Xy. These terms describe the dy- 
namics of an anharmonic oscillator of characteristic fre- 
quency = Y^^^77^/(27r). The next three terms corre- 
spond to the Hamiltonian of a fictitious particle of mass 
^± = ^oC± which is trapped in an anharmonic potential 
along the orthogonal direction X± . These terms describe 
the dynamics of an anharmonic oscillator of characteris- 
tic frequency u± = y^k±/m±/{27T). The anharmonicity 
of the two oscillators is due to the cubic term which re- 
sults from the non-linearity of the Josephson junction. 
The three last terms of the dc-SQUID Hamiltonian mix 
the operators of the two oscillators and consequently cou- 
ple them. Finally, the two-dimensional dynamics of the 
SQUID is similar to the dynamics of two coupled, one- 
dimensional oscillators. We proceed by introducing the 



following dimensionless operators: Xy = y hvp/2E^^X\\^ 



2EyhiypPii, n = Jhv^/2E^Y^ and Px 



/ 2E^ I hu± P± , which verify the commutation rela- 
tions [X||,P||] = i and [X^,Pi] = i. With these new 
operators, the dc-SQUID Hamiltonian can be rewritten 



Hs 



Xn-ahi^„X^ 



+hv:YlX\\ + hvlY^X( + hvlP^P\\ , (18) 

where the parameters a and correspond to the relative 
amplitude of the cubic term compared to the quadratic 
term and hence are direct measure of the degree of an- 
harmonicity of the oscillators. The energies hv^^ hu^ and 
hu^ are the coupling energies between the two oscillators. 

Hereafter, we suppose that the particle is trapped in 
one of the wells of the family [0]. The geometry of this 
well varies as a function of the bias point (Ib^^s) of 
the circuit but does not change with the gate voltage Vg 
and the fiux Therefore, the different parameters of 
the dc-SQUID Hamiltonian only depend on I^j and ^s- 
Fig. [5] shows numerical calculations of the parameters 
of the Hamiltonian under various biasing conditions, us- 
ing the experimental parameters of the circuit studied in 
Ref. [1] (see Appendix [B|) . Figs, [sja) andjsjb) show the 
dependence of the transverse (z/p) and orthogonal 
frequencies as a function of the bias point. Generally, 
the frequency u± is always higher than 35 GHz and is 
at least twice higher than Up. The frequency Up tends 
towards zero when the bias point approaches the critical 
current line. Figs. [5|c,d,e) show the dependence of the 
parameters of the Hamiltonian as a function of the fiux 
for a fixed bias current of 1.89 /iA. Fig. jsjd) shows 
the two anharmonicity parameters of the two oscillators. 
The anharmonicity a is typically around 3 % and in- 
creases close to the critical current line. The parameter 
(j± is very small regardless of the bias, which leads to us 
to the conclusion that the orthogonal oscillator can be 
considered as a harmonic one. Fig. |5|e) shows the differ- 
ent coupling energies. The coupling frequency is of the 
order of 10 MHz and depends only weakly on <l>5'. The 
coupling frequencies and are generally much higher 
than and depend on ^5. Note that vanishes close 
to $5 = and is always negative. Numerical values 
of the parameters of the Hamiltonian, for = 0-1 ^0 
and lb = 1.89 /iA, are given in Tab. [Hj The energy AU 
of the potential barrier (Fig. [5|f)) corresponds to the en- 
ergy which separates the local minimum of the well from 
its closest saddle point. This energy is not a parameter of 
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FIG. 5. Numerical calculations of the parameters of the dc-SQUID Hamiltonian (17), using the experimental parameters of 
the circuit studied in Ref. 1 . The fictitious particle associated to the SQUID is trapped in a potential well of the family [0]. 
(a) Frequencies of the transverse and (b) orthogonal oscillators, i.e. Up and z/^, respectively, as a function of the flux and 
the bias current The black line shows the measured switching current of the dc-SQUID from Ref. [25] which is close to the 
critical current, (c) Frequencies Up and iy± as a function of $s for Ih — 2 /jA, i.e. for the bias points located on the dashed 
line of the Figs, (a) and (b). Measured frequencies us from Ref. [25] are shown as red points, (d) Anharmonicity parameters 
a and a± of the longitudinal and orthogonal oscillators respectively as a function of for = 2 /xA. (e) Different coupling 
frequencies zy*, i/^ and between the two oscillators associated with the terms in y^Xy, y±X^\ and P±P\\ of Hs |l7| . (f) 
Energy AU for the potential barrier as a function of (j)s for = 2 /xA. 



longitudinal oscillator 


transverse oscillator 


coupling 


Up a 
16.24 GHz 3.4 % 


38.26 GHz -0.006 % 


K K K 

-612 MHz -348 MHz 12 MHz 



TABLE II. Numerical values of the parameters of the dc-SQUID Hamiltonian at the working point — 0.1 $o and Ih — 
1.89 11 calculated using the experimental parameters of the circuit studied in Ref. [T]. 



the dc-SQUID Hamiltonian (18). Nevertheless, if A/7 is 



sufficiently small, typically on the order of the expres- 
sion (18) of Hs is too simplified. This occurs when the 



bias point is close to the critical current line. One should 
then take into account the coupling of the quantum levels 



inside the well to those outside the well. Note that this 
coupling is responsible of the escape of the particle from 
the well [26l [27] . This coupling will be neglected in the 
following, assuming that the particle is always trapped 
in a sufficiently deep well. 
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As i^ii , the quantum dynamics of the transverse 

oscihator is much faster than that of the longitudinal os- 
cillator. We will assume in the following that the trans- 
verse oscillator is always in its ground state. It allows 
us to replace in Hs (18) the operators of the transverse 
oscillator by their average values given by {Y±) = 0, 
{YD = 1/2 and {P±) = 0. Only one of the three cou- 
pling terms remains after this simplification. The cou- 
pling reads hu^X\\/2 and can be seen as a modification 
of the bias current of less than 0.5 nA. This term will 
be neglected in the following. Under this condition, the 
dynamics of the particle along the longitudinal direction 
is given by the Hamiltonian 



Vs(X„) 



■Xi 



ahUnXu 



(19) 



In the following, the SQUID dynamics will be studied, 
using this simplified Hamiltonian [26 . We denote by \n) 
and En the eigenstates and the associated eigenenergies 
of Hs^ respectively, such that Hs\n) = En\n), where 
n is an integer number larger or equal to zero. If the 
anharmonicity is weak (a <C 1), the energies En are 
given by a straightforward perturbative calculation and_ 
we find En = {n^l/2)hup-15/4a'^{n^l/2)'^hup. Fi; 



shows the approximate potential of the dc-SQUID 
the three first eigenenergies. When the anharmonicity of 
the dc-SQUID is sufficiently large, the dynamics of the 
dc-SQUID in the presence of an external microwave per- 
turbation involves only the two first levels |0) and |1) [28 . 
In that case, the dc-SQUID behaves as a qubit. Since for 
the dc-SQUID the Josephson energy is much larger than 
the charging energy, the fiuct nations of the phase X\\ are 

much smaller than those of the charge Py. For this rea- 
son, the dc-SQUID is referred to as a phase qubit [14 . 
The dc-SQUID Hamiltonian can be rewritten in the basis 
(|1), |0)), using the Pauli matrices (See Appendix [D|), as 
Hs = hvs/2(7f where vs = (^i —Eo)/h is the character- 
istic qubit frequency. The frequency us can be approxi- 
mated in first order with respect to the anharmonicity a 
as hus = hvp{l — (15/2)cr^). We see that the frequency 
TJs equals the plasma frequency Vp if the dc-SQUID an- 
harmonicity is zero and decreases with increasing anhar- 
monicity (see Fig. jsjc)). When the anharmonicity a is 
zero, the dc-SQUID behaves as a harmonic oscillator de- 
scribed by the typical Hamiltonian Hs = hvp{aa) + 1/2), 
where a) and a are the one-plasmon creation and anni- 
hilation operators, respectively. 



IV. ASYMMETRIC COOPER PAIR 
TRANSISTOR (ACPT) 

The dc-SQUID Hamiltonian is only one of the terms 
that appear in the full Hamiltonian of the circuit. It is 
coupled to a second term associated with the asymmet- 
ric Cooper Pair transistor (ACTP). This section is dedi- 
cated to the theoretical analysis of the ACPT dynamics. 




X„ 



FIG. 6. Approximate one-dimensional potential of the dc- 
SQUID in the direction X\\ . Eigenenergies of the first quan- 
tum states |0), |1) and |2). 



neglecting its coupling with the dc-SQUID. We will first 
build and analyze the ACPT Hamiltonian and will show 
that the ACPT can be viewed as a charge qubit. We 
will then describe the ACPT by using only two charge 
states. The errors on the level of the eigenenergies and 
the eigenstates induced by this simplified description will 
be estimated. 



A. ACPT Hamiltonian 

The ACPT Hamiltonian is identified by isolating from 
the full Hamiltonian of the circuit (15) the terms which 



only contain the operators and n. After some straight- 
forward algebra, the ACPT Hamiltonian reads 

Hacpt =E^ (n - Ugf - EJ cos{S/2) cos{tp - 8/2) 

-^11 sin((5/2) sin(Vi - 5/2)), (20) 

where E^ = E^^ + E^^ and E^ = (2e)V(2CT) are 
the Josephson and charging energies of the ACPT, re- 
spectively. It is the generalization of the Quantronium 
Hamiltonian [16^ for which asymmetries in critical current 
and capacitance of the Cooper pair transistor were ne- 
glected. The gate-charge Ug = CgVg/(2e) corresponds 
to the number of Cooper pairs induced by the volt- 
age applied to the gate-capacitance. The charge states 
\n2e) are the eigenvectors of the charge operator n, such 
^|^2e) = ^|^2e) whcrc Ti = 0, ±1, ±2... is the number of 
excess Cooper pairs on the transistor island. Using the 
commutation relation [V^, n] = i, we identify the action of 

the operator e*'^ = |^2e)((^ + l)2e| which decreases 
the number of Cooper pairs on the island by one unit. In 
the charge representation, the ACPT Hamiltonian can 
be written as 



Hacpt 



--El{n-ng)^ 



^e-(^/^+^)|n2e)((n- 



1 



-^e^(^/2+^)|(n + l)2e)(n2e| 



with £^J/2(cos(V2) +i/isin(V2)) = pj{5)/2e'^, tanx 



2e| 



(21) 



/itan((5/2) and pj{S)'^ 



(cos2(V2) + fi^sm^{6/2)). 
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FIG. 7. Energy spectrum of the ACPT as a unction of Ug 
in the case of a zero Josephson couphng (dashed hues) and 
in the case of pj{6 = 0) = 21.08 GHz (full lines). These 
energy spectra have been calculated numerically, using the 
parameters from Ref. [T]: EJ = 21.08 GHz, = 26.76 GHz 
and /i = —41.6%. The states |— ) and |+) associated to the 
two lowest energy bands correspond to the states of the charge 
qubit. 



The ACPT Hamiltonian is composed of a charging and a 
Josephson term which are proportional to and pj {6) , 
respectively. Let us focus first on the case of a zero 
Josephson coupling {pj = 0). The eigenstates of the 
ACPT Hamiltonian are then the charge states \n2e) with 
the associated eigenenergies E^{n — Ug^. Fig. m shows 
the energy spectrum of the ACPT as a function of Ug for 
a charging energy Eq = 1.28 K. This spectrum consists 
of a series of parabolas, each parabola being associated 
with a specific charge state \n2e) with a minimum en- 



ergy for Ug 



Notice that when the energy difference 



between the ground charge state and the first excited 
charge state is much larger than /c^T, the charge on the 
island is well quantized leading to the Coulomb block- 
ade phenomena [29]. For Ug = 0.5, the energy parabolas 
of the states |02e) and \l2e) cross each other and the 
states |02e) and \l2e) are degenerate. This degeneracy 
is lifted by the Josephson term which couples the neigh- 
boring charge states to each other. The amplitude of 
this Josephson coupling is given by pj . Fig. [S] shows the 
dependence of pj on S for three different Josephson asym- 
metries. Since pj is 27r periodic in 6, the range of 6 has 
been restricted to the interval between — 27r and 27r. In 
the case of a symmetric transistor {p = 0), the Joseph- 
son coupling is maximum for ^ = 0, equal to Ej; it is 
zero at (5 = ±7r. For a finite asymmetry /i, the Josephson 
coupling reaches a maximum for J = and equals Ej^ 
and a minimum equals to pEj for S = ±7r. For a Cooper 
pair box {p = ±100%), the Josephson coupling does not 
depend on S and remains equal to Ej. The Josephson 
coupling then depends strongly on the Josephson asym- 
metry, especially for S = ±7r where it can vary from zero 
to EJ. The full energy spectrum of the ACPT, which 
takes into account the Josephson coupling, is calculated 
numerically by diagonalizing the Hamiltonian (21) with 
8 charge states. It is plotted as a function of Ug in Fig. [7| 



FIG. 8. Josephson coupling pj as a function of the su- 
perconducting phase 6 across the ACPT, for three different 
Josephson asymmeties: p — for a symmetric transistor, 
p = —41.5% for an asymmetric transistor and p = 100% for 
a Cooper pair box. 



for a fixed superconducting phase ^ = and for the 
parameters of [J: = l.2^K, Ej = IMhsK and 
p = —41.6%. The spectrum consists of several energy 
bands which do not cross each other, leading to an en- 
ergy gap between the lowest and the first bands. We 
associate with these two bands the ground state | — ) and 
the first excited state |+), respectively. These two states 
correspond to the states of a qubit, with a characteristic 
energy hur- Note that to operate the ACPT as a qubit, 
the gate-charge should be close to 0.5 (modulo 1) where 
the transition energy from the state |+) to the third level 
is higher than the frequency ut- Indeed, when the ACPT 
Hamiltonian (21) is subject to an adequate perturbation. 



the quantum dynamics of the ACPT will only involve the 
states |— ) and |+). Because Ej > E^ the states |— ) and 
1+) can be written as a superposition of a few charge 
states (typically four). For this reason, the qubit formed 
by the ACPT is referred to as a charge qubit. The Hamil- 
tonian of the charge qubit can be written using the Pauli 
matrix (see appendix IdI) as 



2 ^ 



(22) 



The frequency vt of the charge qubit depends on the 
two parameters 8 and Ug as shown in Fig.joj^a). This fre- 
quency, as is Pj, is also 27r periodic as a function of 5. It is 
maximum and minimum at the points {5 = 0,n^ = 1/2) 
and {5 = 7r,n^ = 1/2), respectively. These are optimal 
points where the ACPT is, in first order, insensitive to 
charge, flux and current fluctuations. Fig. [lO] shows the 
experimental frequency of the charge qubit measured in 
Ref. [1] as a function of 6. The theoretical fit, shown 
in red, is obtained by diagonalizing the ACPT Hamil- 
tonian (21) in a basis of 8 charge states. It allows us 



to accurately find the parameters of the charge qubit, 
i.e. the Josephson energy Ej = 21.08 GHz, the charging 
energy E^ = 26.76 GHz and the Josephson asymmetry 
p = —41.6 %. Nevertheless, the capacitance asymmetry 
can not be extracted from the fit since it does not enter in 
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FIG. 10. (Color online). Frequency i^t as a function of the 
phase 6. The points are the experimental data from Ref. [l]. 
The solid curve is the theoretical frequency ut calculated by 
diagonalizing the ACPT Hamiltonian (21) with the parame- 
ters Ej = 21.08 GHz, = 26.76 GHz and /x = -41.6 %. 
Inset: Frequency ut as a function of Ug for S = 0.5497r. The 
points are the experimental frequencies ut and the solid line 
is the theoretical frequency ut calculated as before. 



function of the charge states as 



'a) Frequency ut of the charge qubit 



FIG. 9. (Color online), 
as a function of Ug and 6 for the experimental parameters of 
Ref. P: E^ = 1.28 /c^K, EJ = 1.01 ksK and /x 



(b) Frequency ut as a function of rig for ^ = 0, 

S = 7T. 



= -41.6%. 

S — ±7r/2 and 



the ACPT Hamiltonian (21). The Josephson energy and 



the capacitance of a Josephson junction of the ACPT are 
both proportional to the junction surface and, therefore, 
the capacitance and Josephson asymmetries are equal in 
first approximation. We will see in Sec. [V| that the ca- 
pacitance asymmetry can be extracted from the coupling 
between the dc-SQUID and the ACPT and we will find 
A = 0.875/i. 



B. Description of the ACPT with two charge levels 

If the Josephson coupling pj is much smaller than the 
charging energy and if rig ^ 0.5, the qubit states 
can be expressed as a superposition of the two charge 
states |02e) and |l2e)- In that case, the Hamiltonian of 
the transistor is simply given by its matrix form H^^pj. 
which reads, in the charge basis (|02e),|l2e)), 



ACPT 



-pjC 

W2+x)/2 El{l 



i(S/2+x) 12 
2 



n„ 



(23) 



The 



Hamiltonian 

(l-n,)2) ± 

2 i/^5^(l — 2n^)2 + p^, and the qubit energy reads 



eigenvalues 
are given 



by El 



of this simplified 



El 



E^_ 



Note that for n„ 



0.5, we have 



Pj. The eigenstates |— o) and |+o), associated 
with the energies E- and E-^., can be written as a 



I+O) =«*|02e)+/3|l2e) 
I-O) =-/?*|02e) + a|l2e) 

7r)/2 



(24) 



2|pjl 



with a = cos(6'/2)e'(''/2+x- 
sin(6'/2)e'('^/2+x-7r)/2 ^^^g 

(^e [0,7r]). 

For the ACPT studied experimentally in Ref. [1 , 
condition Pj{S) <C E^ is too strong, especially at S 
where the ratio pj/E^ ^ 79% is maximum. For this rea- 
son, we would like to quantify the error induced by the 
description of the ACPT with only two charge states. Let 
us first focus on the error in the qubit energy. Fig. pT 



-2nJ 

the 
= 



shows the difference between the frequencies vj. calcu- 
lated with two charge states and the real frequency vt as 
a function of (5 for = 0.5. This difference is minimum 
for 5 = ±7r and equals 58 MHz, corresponding to an error 
in energy of 0.5%. For (5 = 0, the difference is maximum 
and equals 780 MHz, corresponding to an error in energy 
of 3.8%. 

We now discuss the error made on the level of 
the states |— o) and |+o)- The non-reduced Hamil- 
tonian of the ACPT can be^rewritten ^s a function 
of Hacpt Hacpt = E[acpt + where W 
perturbative term which takes the form W = 

ny^{0,l} 

theory, we calculate the probability |(— o| — )P (l(+o|+)P) 
that the state |— ) (|+)) is in the the state |— o) (|+o))- 



IS a 

2 



-ng)'^\n){n\. Using first-order perturbation 



Fig. [TT] shows the dependence of these probabilities as a 
function of S for rig = 0.5. The probability |(— o|— )P 
is always smaller than the probability |(+o|+)P- In- 
deed, the energy of the state |+o) is closer to the charg- 
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800 





FIG. 



12. (Color online). 

|2 



Probabilities |(- 



-)|^ (black) and 



|(+o|+)|^ (blue) as a function of the ratio / for Ug — 0.5 
and ^ = 0. Numerical calculations of these probabilities (full 
lines) are compared with first-order perturbation calculations 
(dashed lines). The vertical dashed line indicates the values of 
these probabilities at the ratio E^ / E^ = 79% corresponding 
to those of the ACPT studied in Ref. [T]. 



COUPLING 



FIG. 11. (Color online). Errors in the energies and states of 
the charge qubit of Ref. [1 induced by describing it with only 
two charge states, (a) Difference between the qubit frequency 
calculated with two charge states and its real frequency 
vt as fonction of 5 for rig — 0.5. (b) Probabilities |(— o|— )|^ 
(black) and |(+o|+)|^ (blue) as a function of 5 for rig — 0.5. 



ing energies of the states | — l2e) and |22e) and conse- 
quently more perturbed by these two states. The prob- 
abilities are minimum for 5 = (|(— o|— )P = 97.4% and 
|(~o|~)P = 94.7%) and, therefore, for this value of the 
phase, the error made by considering the states of the 
ACPT as |— o) and |+o) is maximum. Fig. pTj shows the 
dependence of the probabilities |(— o| — )P and |(+o|+)P 
as a function of the ratio Ej/E^ for ^ = 0. Theses 
probabilities have been calculated with a first-order per- 
turbation theory and also numerically by using 8 charge 
states. We see that when Ej/E^ < 1, the first-order 
perturbation theory agrees quite well with the numer- 
ical simulations. But when Ej/E^ > 1, the analytical 
calculation gives probabilities significantly lower than the 
numerical one. The probabilities |(— o|— )P and |(+o|+)P 
decrease with increasing Ej /E^, which is explained by 
the fact that the Josephson coupling mixes more and 
more the charge states |02e) and |l2e) with the other, 
closer-in-energy charge states. In order to simplify the 
calculation of the analytical expression of the coupling 
(see below), we consider in the following that | — ) = |— o) 
and 1+) = 1+0)7 but do not approximate the qubit fre- 
quency to u^. 



So far, we have considered independently the quantum 
dynamics of the longitudinal mode of the dc-SQUID and 
the ACPT. However, in the studied circuit the dc-SQUID 
and the ACPT are connected in parallel and therefore 
coupled to each other. The independent dynamics of 
the ACPT and the dc-SQUID has to be reconsidered 
especially when the two qubits are close to resonance 
{i^s ^ i^t)- In this case the coupling effects are the 
strongest. In this section, we derive the expression of 
the coupling Hamiltonian by considering the ACPT as a 
charge qubit and the dc-SQUID as either a tunable har- 
monic oscillator or a phase qubit. We will see that the 
total coupling is the sum of two distinct contributions: a 
capacitive and an inductive Josephson coupling. 



A. Capacitive coupling 

The capacitive coupling Hamiltonian couples by defini- 
tion the charge of the dc-SQUID with that of the ACPT. 
It reads 



Hcoupl,Capa = -^^1/ ^T^lK^ ~ ^s)' (^5) 

^n|| y 2E^^ 

We consider hereafter two different limits for the dc- 
SQUID in order to simplify this capacitive coupling 
Hamiltonian. 

The first limit corresponds to a dc-SQUID with an 
anharmonicity factor a equal to zero. This limit can 
be achieved when the fictitious particle associated with 
the dc-SQUID is trapped in a deep well, which is gen- 
erally true when the dc-SQUID is biased at the working 
point lb = and = 0. Under this condition, the 
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dc-SQUID behaves as a harmonic oscihator and several 
levels are involved in the dynamics. The momentum op- 
erator in the charge coupling Hamiltonian is then given 
by P|| = {d^ +^)/ a/2, where d^ and d are the one-plasmon 
creation and annihilation operators. If we describe the 
charge qubit with two charge states (Eq. (|24])), the opera- 
tor h can be written as n = ^ /+ ^ sin(^)(jj — ^ cos(^)aJ, 
where / is the identity operator and and are 
the Pauli matrices defined in the eigenstate basis of the 
charge qubit. For Ug = 1/2, we have 6 = 7r/2 and the 
charge coupling Hamiltonian simplifies to 



Coupl,Capa 



(26) 



where Ec^c = ^ /Cn\\\j hvp/ E^^ determines the strength 
of the capacitive coupling. 

The second limit is realized when the anhamonic- 
ity cr is typically around 3 % or larger; then the dc- 
SQUID can be described as a phase qubit. Generally, 
(J does not exceed 10% (see Fig. |5|d)) so the two low- 
est eigenstates of the dc-SQUID Hamiltonian are very 
close to the two lowest eigenstates of the harmonic os- 
cillator. Consequently, the momentum operator can be 
expressed in terms of the eigenbasis of the phase qubit as 
P|| ^ (|0)(1| + |l)(0|)/\/2 = af/\/2 . For Ug = 1/2, the 
charge coupling Hamiltonian takes the following form: 



Hcoupl,Capa = -Ec^cO-^O- 



(27) 



The capacitive coupling produces a transverse interac- 
tion between the two quantum systems. This coupling 
was already discussed in Ref. [30, 31 where a Cooper 
pair box was coupled to a harmonic oscillator. The ca- 
pacitive coupling is vanishing in the limit of small Cj. 
We notice also that £^c,c depends on the bias variables 
through and therefore is weakly tunable. 



basis: sin('^ 
ing Eq. (|24| 

comes 



-S) = ^ (e^^|l2e)(02e| - e"^' |02e) (l2e |) • Us- 

the Josephson coupling Hamiltonian be- 



Coupl,Jos 



where Er 



-—iEcj{d - a"^) [sin((5/2 - x)(sin((9) al 
+ cos{0) aj) — cos{5 /2 — x)o; 



.T1 

y J ' 



(29) 



= a\\/2Ej2yE^^/hvp quantifies the 
strength of the Josephson coupling. For rig = 0.5, we 
have 6> = 7r/2 and the Josephson coupling Hamiltonian 
reduces to 



Hcoupi.Jos =iEcj{d - d^){cos{S/2 - x)^y 
-sin(V2-x)^J). 



(30) 



In the limit of finite anharmonicity cr, the dc-SQUID 
can be approximated by a two-level system. In that case, 
the position operator can be written in the eigenbasis of 
the phase qubit as X\\ ^ i(|0)(l| - |1)(0|)/V2 = cr^/y/2. 
Using the latter expression, the Josephson coupling be- 
tween the charge and phase qubits takes the following 
form for r 



1/2: 



Hcoupl.Jos =Ec,j cos((5/2 - x)cry cr. 



-Ec,jSm{6/2-x)^y^, 



(31) 



The Josephson coupling contains two different terms. 
One describes a transverse coupling cTy&y or (a — d^)ay 
which gives rise to coherent energy exchange at the reso- 
nance between the charge qubit and the phase qubit (or 
the oscillator). The effects of this first term on the quan- 
tum dynamics of the circuit are similar to those produced 
by the capacitive coupling term in or {d — d^)a^ . 

The second term (J^aJ or (a — d^)a^ contains a trans- 
verse contribution for the SQUID and a longitudinal term 
which depends on the transistor qubit state. Its contribu- 
tion will explain the quantum measurement of the charge 
qubit in the us <C i^t limit (see Sec. |VI C" ). We notice fi- 
nally that the two terms are strongly tunable with the 
bias parameter S. 



B. Josephson coupling 



The Josephson coupling Hamiltonian couples the 
phases Xy and lip associated with the dc-SQUID and the 
transistor, respectively. This coupling is mediated via 
the Josephson junction 2 of the ACPT and reads 



Coupl,Jos 



12E, 



c 



hi 



X\\sm{iP-S). (28) 



We derive now the expression of the Josephson cou- 
pling Hamiltonian for the two different limits of the dc- 
SQUID. When the dc-SQUID is in the harmonic oscilla- 
tor limit, the position operator can be written as Xy = 
i{d — d^)/\^. On the other hand, close to = 1/2, the 
operator sm{ip — S) takes the following form in the charge 



VI. QUANTUM DYNAMICS OF THE 
COUPLED CIRCUIT 

The full Hamiltonian of the coupled circuit is given 
by the sum of the Hamiltonians of the dc-SQUID, the 
ACPT and the coupling. It reads: H = Hs + Hacpt + 
Hcoupi,Capa-\-Hcoupi,Jos- In Order to simplify this Hamil- 
tonian, we consider hereafter the situation where the 
gate-charge is fixed to Ug = 0.5. This charge value has 
been mainly used for the charge qubits experiments (see 
Refs. [H HSl [32] ) . Indeed, at Ug = 0.5, the charge qubit 
is insensitive in first order to charge noise and, therefore, 
its coherence time is longer. We will first derive the ex- 
pression of the full Hamiltonian when the two quantum 
systems are close to the resonance condition. We will 
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consider the dc- SQUID either as a phase qubit or as a 
harmonic oscihator. Finahy, we wih discuss the quantum 
measurement scenario of the transistor states. Its princi- 
ple derives from the couphng to the SQUID, enabUng to 
use it as a detector of the charge qubit 's state. 



A. A charge qubit coupled to a phase qubit 

We first consider the dc- SQUID as a phase qubit such 
that the full Hamiltonian governs the quantum dynam- 
ics of two-coupled qubits. Using the expressions of the 
Josephson and capacitive coupling in that limit, we find 



H 



'S .S 



S -Tf 



cos((5/2 - x)^y^i - Ec,j sin(J/2 - x)^y^i (32) 

Let us introduce the raising operators and and the 
lowering operators and a^. These operators are de- 
fined by (j^ = ((jj ± iay)/2 and = (af ± icry)/2. 
The terms <3-faJ and cfyCfy can be written using the 
four products cf^^^^cf^cf^ and a^a^^ whereas the 

term <J^aJ is a function of the operators d-faj and d-^aj . 

The product afa^ corresponds to an excitation of the 
phase qubit and de-excitation of the charge qubit. This 
coupling term mediates the transition between the states 
|0, +) and ). This transition is only relevant close 
to resonance, where it contributes to the low-frequency 
dynamics of the coupled system. Far away from reso- 
nance, this term gives rise to high-frequency dynamics 
(frequencies of the order of vs and ut) which is aver- 
aged away on the typical time scales of the experiment 
[1 . The rotating- wave approximation consists in neglect- 
ing these non-resonant terms, called also inelastic terms. 
The previous reasoning applied to the term leads 
to the same result. The products d-fa^ and a^a^ couple 
the states |0, — ) and |1, +). The transition between these 
two states leads to high-frequency dynamics and are ne- 
glected hereafter. For the same reason, the coupling term 
o-yO-J will be ignored. Finally, close to the resonance, the 
full Hamiltonian simplifies to 



H 



'S .s 



S I ±.T\ 



(33) 



where the coupling strength g is given by ^ = 2£^c,c — 
2Ecj cos{6/2 — x)' We assume hereafter that the gate 
capacitance Cg <C and <^ Co which is true for 
the measured circuits in Refs. [l]|16]. In these limits, we 
have 1 1 ^ 2Co(C2^+Cf )/(a||Cj). The coupling energy 
can then be written as a function of the capacitance (A) 
and Josephson (/i) asymmetries of the ACPT as: 



g[(l + A)/..,-(l + 



fi)Ejcos{S/2-x)]- 
(34) 



Further simplification of the coupling energy can be ob- 
tained using the fact that the charge qubit is described 
with two charge states, and that the SQUID has a zero 
anharmonicity. Then, the transistor frequency is related 
to S by = {Ej/h){cos\S/2) ^ fi^sm\S/2))^/^ (see 
Sec. IV B) and we have Up = vs- By using these rela- 
tions, the coupling g (l34|) can be rewritten as 




(l+A)i/5-(l + /i)^T(cos2(x) 



sin^(x)> 



(35) 



Note that coupling g depends strongly on the bias param- 
eters, via the phase 5 and the frequency vs- Therefore, 
the proposed circuit presents an intrinsic tunable cou- 
pling between a charge and a phase qubit, which we will 



analyze in more detail in Sec. VII This transverse cou- 
pling enables t o realiz e two-qubit gates operation as for 
example the V iS WAP operation |33] . 



B. A charge-qubit coupled to a tunable harmonic 
oscillator 



The Hamiltonian (33), which describes the dynamics 



of two coupled qubits, is valid if the second and highest 
levels of the dc- SQUID do not participate in the quantum 
dynamics of the circuit. This is the case when the an- 
harmonicity of the SQUID is sufficiently strong. In the 
limit of zero anharmonicity, the quantum dynamics of 
the coupled circuit is described by the Jaynes-Cummings 
Hamiltonian [34 



H = hup[aa^ ^ 2^ ^ ~Y~ ^ 



T ^(a^a^^a&l). (36) 



This Hamiltonian (36) is very similar to the one ob- 
tained by coupling a charge [35] or a transmon [36 or 
phase qubit [37 to a coplanar waveguide cavity. In our 
circuit the resonator is realized by a micro dc-SQUID. It 
can be viewed as "micro-resonator" more convenient for 
integration than usual coplanar resonators since its size 
ranges in the micrometer scale three orders of magnitude 
smaller. Moreover its resonance frequency is strongly 
tunable. However up to now it suffers from a much 
shorter coherence time. 



C. Quantum measurements of the charge qubit by 
the dc-SQUID 

Our theoretical analysis introduces two different kinds 
of coupling which can be used to perform quantum mea- 
surements of the charge qubit by the dc-SQUID. Here 
we will apply this analysis to the quantronium read-out 
|T^ and the adiabatic quantum transfer method used in 
Ref. p. 
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1. Quantronium read- out 

This read-out is obtained the hmit where vt ^ i^s- As 
we win see, even if the qubits are far from resonance the 
couphng stih affects the dynamics of the circuit. Indeed, 
since the dynamics of the ACPT is much faster than the 
dc-SQUID dynamics, the dc-SQUID is only sensitive to 
the average value of the ACPT operators. After some 
straightforward algebra, the effective Hamiltonian of the 
dc-SQUID takes the following form: 



l±) 

S,eff 



Cr^W dUa 



dEi 



l±> 



dS 



Xu 



(37) 



with E\_^^ and E\_^ the eigenenergies of the ACPT associ- 
ated with the states |+) and |— ). For Ug = orn^ = 1/2, 
we have dE\j^^ / dug = and the SQUID Hamiltonian 
simplifies to 



LU 



— asymmetric CPT 

— symmetric CPT 

— Cooper pair box 




FIG. 13. (Color online). Effective current r~^l at Ug — 1/2 
as a function of the phase 5 for an asymmetric (/i = 0.42) and 
symmetric (/x = 0) transistors and a Cooper pair box (/i = 1). 



l±> 

S,eff 



(38) 



where / 



l±> 

add 



-{a\\/(j)o)dE\±^ /d5 adds to the bias current 
lb and depends on the state of the ACPT | — ) or |+). For 
Ug = 1/2, in the limit of two charge states only, the 
additional current reads 

lining = 1/2) = ±^(1 + m) sin(<5/2 - x)- (39) 



The currents /i+^(ng = 1/2) and ll^^{ng = 1/2) take 
opposite values and, therefore, can be used to determine 



the transistor state. Fig. 13 shows the normalized cur- 
rent iI^J^/{EJ /2(I)q) for Ug = 1/2 as a function of the 
phase S for an asymmetric {ja = 0.42) and a symmet- 
ric {fi = 0) transistor and a Cooper pair box (/i = 1). 
For a symmetric transistor, the current has a strong dis- 
continuity at (5 = ±7r jumping between the two extreme 
values —Ej/2(j)o and Ej /2(j)o. This discontinuity dis- 
appears for a finite Josephson asymmetry; its maximum 
value is lower. Naturally, for a Cooper pair box, the ad- 
ditional current is zero as one of the junctions is replaced 
by a pure capacitance. 

As the switching probability of the dc-SQUID depends 
strongly on the bias current [26 the additional current 
ladd can be used to detect the ACPT state. For instance, 
in the Quantronium circuit [16 , the state of a symmet- 
ric Cooper Pair Transistor (CPT) is read out by mea- 
suring the switching probability of a Josephson junction 
placed in parallel with the CPT and which replaces the 
dc-SQUID in the studied circuit. In that case, the charge 
qubit read-out is explained by the X^aJ term resulting 



from the Josephson coupling (28). 

The expression ( 39 ) of the additional current has been 
established under the condition ut ^ ^s- Is this ex- 
pression still valid when ut ^ ^5? In order to answer 
this question, we compare experimental data and the- 
oretical predictions for the difference of the additional 
currents at Ug = 1/2 and = in the ground state 
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FIG. 14. (Color online), (a) Experimental (black squares) and 



theoretical (red line) current difference A/^^^ = l\^^(ng 



Ildd('^9 = 0) as a function of the flux $5 (or the phase 



1/2)- ^1"^ 

S indicated in the upper scale), (b) Frequencies ratio vt/j^s 
as a function of $5 



of the ACPT. This difference noted A/ 
rl 



1 I 



as a function 



1/2) - lidding = 0) is shown in Fig. 

of the flux ^s- The theoretical curve calculated with- 
out any free parameters agrees well with the experiment 
when vt/^s ^ I7 but differs from the experiment when 
^ ^s- The measured current amplitude drops when 
vt is close to vs suggesting a drop in the contrast of the 
quantum measurement when ut ^ ^s- 
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FIG. 15. (Color online). Spectroscopies of the phase (blue 
curve) and charge (red curve) qubit studied in Ref. 1 at (a) 
Ug — 0.5 and (b) Ug ~ 1 (vt ^ ^s) for the working point 
(/^^^^ = 1957 nA,$^^ = -0.064 $o). These spectroscopies are 
measured with a nanosecond flux pulse which changes the flux 
in the SQUID loop from to (c) Evolution of the 

qubit frequencies vt and vs during this pulse. Blue and red 
points indicate the frequencies vs and vt at the working point 
of spectroscopy (a). Inset shows an antilevel crossing in the 
energy spectrum where an adiabatic transfer happens during 
the flux pulse, the state |1,— ) (|0,+)) being transferred to 
the state |0,+) (|1, -)). 



2. Adiabatic quantum transfer 

In Fig. [Tsj the first energy levels of the charge and 
phase qubit in the coupled circuits are plotted. The quan- 
tum measurement of the charge qubit is performed by a 
nanosecond flux pulse which transfers the quantum state 
|0, +) prepared at the working point (I^^,<l>^^) to the 
measurement point (I^^^,<l>|^^). At that point the SQUID 
is very close to the critical line and the escape probabil- 
ity is finite. Spectroscopy measurements show clearly the 
read-out of the charge qubit by this method even in the 
limit Vt ^ 'i^S' How can we explain this read-out? 

The coupling terms cr^a^ and cr^a^ produce an anti- 
level-crossing the amplitude of which depends on the 
coupling strength g (see Inset of Fig. [l5|c)). During a 
nanosecond flux pulse, the coupled system remains in 
its original energy state; as a result, the quantum state 
|0, +) evolves adiabatically into to the state |1,— ). Due 
to large value of the coupling strength close to the es- 
cape point, Landau-Zener transitions can be neglected. 
The final state, i.e. |1, — ) or |0, — ), is determined by the 
switching measurement. 



VII. TUNABLE COUPLING 

In this section, we calculate numerical values of the 
coupling strength ^, using the circuit parameters of 
Ref. [1 and show that the coupling can be tuned over 
a wide range with the bias parameters. 



FIG. 16. (Color online). Color plot of the analytical cou- 
pling ( |35[ ) as a function of the bias current Ih and the flux 
for the wells family [0] of the dc-SQUID. The sohd hne 
indicates where the two qubits are in resonance {vs — ^t). 




FIG. 17. Coupling g derived from Eq. (35) as a function of 
5 when the two qubits are in resonance [vt — vs) for three 
different charge qubits: (red curve) an asymmetric transistor 
with A = /i = —41.6%, (bue curve) a symmetric transistor 
\ — 11 — and (black curve) a Cooper pair box with A = 
—41.6% and /x = — 1. The calculations have been realized by 
using the parameters of Ref. [1 given in Appendix p| 



a function of the bias current /5 and the flux for the 
family [0] of the dc-SQUID. We find that the coupling 
between the two qubits can be tuned from zero to more 
than 1.2 GHz. This on-off coupling is one of the needed 
requirements to realize ideally one- and two-qubit gates. 
In Ref. [1 , the coupling g has been measured at reso- 
nance, where the coupling effect on the qubits is maxi- 
mal. Far away from resonance, the eigenenergies of each 
individual qubit are shifted by the amount ^^/4| A| where 
A = hvs — hi^T is the detuning. So, if the detuning is 
large the coupling can be neglected. The solid line of 
Fig. [16] shows where the qubits are in resonance and, 
consequently, where the coupling can be easily measured 
by spectroscopy. 
Fig. 



Fig. 16 presents the dependence of the coupling g as 



T7| shows the dependence of the coupling ^ as a 
function of the phase S for three different charge qubits: 
an asymmetric transistor fll, a symmetric transistor [16] 
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and a Cooper pair box \3T . The coupling has been calcu- 
lated using the expression ( 35 ) at the resonance us = 
with the parameters of Ref.^lT and only the ACPT asym- 
metries have been varied. For the asymmetric transis- 
tor with A = /i = —41.6%, the coupling g (red curve) 
is maximum at (5 = ±7r where it equals to 1217 MHz; 
it becomes zero at (5 = 0. In the case of the symmet- 
ric transistor with A = /i = 0, the coupling strength 

reads g = Oi\\/2\J E^jjvs {ys — ^t)- Consequently, at 
the resonance, the coupling between a symmetric tran- 
sistor and a SQUID is zero. For a Cooper pair box with 
(/i = — 1 and A = —41.6%), the coupling reads 



^J2 



g = E^hvsil + A). This corresponds to the result 
of Ref. [3j . The calculated coupling g (black curve) does 
not depend on 5 and remains equal to 514 MHz. 

Fig. [iSj^a) shows the dependence of the analytical and 
numerical couplings at the resonance as a function of 
the absolute value of the phase 5. The numerical sim- 
ulations allow us to check the validity of the analytical 
coupling. These simulations have been realized by diag- 
onalizing the full Hamiltonian in the basis of 20 charge 
states and the first 9 excited states of the dc-SQUID in 
the absence of anharmonicity. The dc-SQUID and ACPT 
Hamiltonians, as well as the Josephson and capacitive 
coupling Hamiltonians can all be expressed naturally in 
these bases. The numerical coupling is found by calcu- 
lating the energy spectrum as a function of 5 between 
— TT and TT for a fixed frequency vs- The spectrum is first 
calculated without the coupling terms in order to find 
the two opposite values oi 5 = ±(5^ where the ACPT and 
the dc-SQUID are in resonance {ys = ^t), i-e. where the 
second and third energy bands intersect. In the presence 
of coupling, the degeneracy between the eigenstates is 
lifted and an anti-level crossing appears with an energy 
separation equal to the coupling g 5 = As an 

example, the inset of Fig. 18 ^a) shows the energy spec- 
trum for (5^ = 0.57r with and without coupling from which 
we extract g/h = 251 MHz. The numerical and analyti- 
cal simulations remarkably give quite close results which 
confirms the validity of the analytical expression (35). 



The theoretical coupling is calculated here without any 
free parameters by assuming equal Josephson and capac- 
itive asymmetries {X = ja = —41.6%). Finally, we note 
the good agreement with the experimental coupling mea- 
sured in Ref. [1 and shown in black points in Fig. 18 By 
adjusting the capacitive asymmetry to A = —36.4% the 
numerical coupling is in perfect agreement with the ex- 
periment as shown in Fig. ITsFb). 



VIII. CONCLUSION 

In conclusion, we have analyzed in detail the theoreti- 
cal quantum dynamics of a superconducting circuit based 
on a dc-SQUID in parallel to an Asymmetric Cooper Pair 
Transistor (ACPT). The Lagrangian of the circuit was 
first established from the current conservation equations 
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FIG. 18. (Color online), (a) Coupling ^ as a function of 5 
when the SQUID and the ACPT are in resonance (ut = J^s)- 
The black points are the experimental couplings measured in 
Ref. [1 . The blue line is the theoretical coupling calculated 
from the analytical expression (35), using the circuit parame- 



ters of Ref. U] and equal Josephson and capacitance asymme- 
tries fi — \ — —46.6 %. The red line is the numerical coupling 
for the same circuit parameters. It is calculated by diagonal- 
izing the full Hamiltonian in the basis of 8 charge states and 
the first 9 harmonic dc-SQUID states. The inset shows the 
numerical simulations of the eigenenergies of the circuit with 
(full curve) and without (dashed curve) coupling from which 
we extract the coupling gjh — 251 MHz at ^ = 0.57r. (b) Nu- 
merical coupling ^ as a function of 8 when the SQUID and the 
ACPT are in resonance, calculated using the circuit parame- 
ters used in (a) but with a different capacitance asymmetry 
A = 0.87/i = -36.4%. 



expressed at each node of the circuit. The Hamiltonian, 
deduced from the Lagrangian, is decomposed in three dis- 
tinct terms, namely the dc-SQUID, the ACPT and the 
coupling Hamiltonians. We first studied the individual 
dynamics of the dc-SQUID and the ACPT. Depending 
on its anharmonicity, the dc-SQUID can be seen either 
as a harmonic oscillator or as a phase qubit, whereas the 
ACPT behaves as a charge qubit. In addition to the opti- 
mal bias points ((5 = 0,n^ = 1/2) which was successfully 
demonstrated in a symmetric Cooper pair transistor, the 
ACPT presents a second optimal point ((5 = tt, = 1/2). 



At these points, the charge qubit is insensitive in first 
order to the charge, flux and current noise and there- 
fore shows a larger coherence time. We found that the 
coupling Hamiltonian between the dc-SQUID and the 
ACPT is made of two different terms corresponding to 
the Josephson and the capacitive couplings, which mix 
phases and charges of both sub-circuits, respectively. 

The Hamiltonian of the full circuit was discussed 
through two different limits of the dc-SQUID. When it 
behaves as a harmonic oscillator, the quantum dynam- 
ics of the circuit is described by the well-known Jaynes- 
Cummings Hamiltonian. The microscopic circuit is then 
similar to a two-level atom coupled to a single-mode op- 
tical cavity. It offers compared to this latter a better 
tunability, a faster control and read-out of the quantum 
system, and a good scalability for complex achitectures 
implementation. For example a circuit of several ACPT 
qubit in parallel could be considered whose quantum in- 
formation will be mediated by the dc-SQUID [39 . When 
the anharmonicity of the dc-SQUID is strong, it behaves 
as a phase qubit. The full circuit is then described as the 
coupling of two different class of qubits, i.e., a phase and 
a charge qubit. The coupling Hamiltonian contains terms 
in d-y&y and o-Jaf which are prominent when the two 
qubits are in resonance. These terms allow two-qubits 
gate operations as the V iS WAP gate. In addition they 
enable to read out the quantum state of the ACPT by a 
nanosecond flux pulse as observed in Ref . [1] . Indeed such 
a pulse produces an adiabatic quantum transfer of the 
state |0, +) into the state |1, — ), i.e. the energy quantum 
stored in the ACPT is transferred into the dc-SQUID in 
order to be detected. A non-resonant term in (a — d^)d-J 
or in o-yO-^ is present in the Josephson coupling. Al- 
though its effect on the quantum dynamics is weak, this 
latter term explains the charge qubit read-out in the limit 
^ by mean of an effective additional current in 
the dc-SQUID. That read-out method is employed in the 
Quantronium circuit [16. . 

In both limits of the dc-SQUID, we demonstrated that 
the coupling can be strongly tuned, mainly with the 
Josephson coupling term which has a strong phase S de- 
pendence. It can be used to accomplish two-qubit gate 
operations, and can also be turned off in order to per- 
form one-qubit gate operations without desturbing the 
unaddressed qubit. 
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Appendix A: Current conservation laws 

The current conservation law, applied to each node of 
the circuit (see Fig. [ijb)), yields six equations for the ac- 
tive phases (pi^ (p2^ and ^. These equations are iden- 
tical to the Euler-Lagrange equations derived from the 
circuit Lagrangian, ([2|, ([3|, Q and read, respectively. 



00 ^ j-^^ - Co(l)o(pi - lo sm{(pi) = 0; 
^1 

(j)o — ^000^2 - ^0 sin (p2 = 0; 



(Al) 
(A2) 



-CgMi^ - 0) - C^M - sin(V^) = 0; (A3) 



(A4) 



7-^ 



-C^Mi - ^) - II sin(^ - V; - (/.x) = 0. (A5) 



In Sec. |II C[ we show that this system of six equations 
can be reduced to four equations by ignoring the high 
frequency quantum dynamics of the phases 7 and ^. 
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TABLE in. Electric parameters of the coupled circuit studied 
by Pay et al. 1 . 



—hn 
-huQ 



dC 



dC 

dC 

dC 

09 ' 



= 4 {(2Co + ajC^Xii - a\\Cl (aj^^ + ^ 
= 4 {(2Co + alCl)Y^ - a^Cl (a\\X\\ - ^ 

4>l[{Cp + Cg)e-Cg^}. 



Appendix B: Parameters of the coupled circuit 
studied in Ref. [T] 

Throughout this article, we illustrate the theory with 
numerical values and plots calculated by using the pa- 
rameters of the circuit studied in Ref. [1 . These param- 
eters are collected in Tab. IIIII 



Appendix D: Pauli matrixes 

The Pauli matrices related to the dc-SQUID are de- 
fined in the eigenbasis of the phase qubit {|1), |0)} as 

^^ = [0-1)^ ^- = [10)^ ''y = [i )■ 

(Dl) 

Similarly the Pauli matrices related to the ACPT (<jj, 
crj, ay) are defined in the eigenbasis of the charge qubit 



Appendix C: Conjugate variables 



The phases (Xy , F^, ^) and their conjugate mo- 
menta {—hP\\^—hP±^—hn^—hnQ) are the appropriate 
variables of the circuit Hamiltonian (|T5|). The momenta 
are related to the velocities(X|| , 1^, ^) involved in the 
kinetic part of the Lagrangian ([l3| by the following ex- 
pressions ^22j : 
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